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^The  asymptotic  stress  and  deformation  fields  of  a  crack  propagating  steadily  and 
quasi-statically  into  an  elastic-plastic  material  are  presented.  The  material  is  characterised  by 
J2  deformation  theory,  suitably  modified  to  account  for  unloading  and  reloading,  together 

with  linear  strain-hardening.  The  cases  of  plane-strain  and  plane-stress  mode  I  are 
considered.  The  associated  governing  equations  are  integrated  analytically,  and  the 
boundary  conditions  lead  to  algebraic  nonlinear  equations  which  are  solved  numerically. 
Results  are  given  for  the  strength  of  the  singularity,  and  for  the  distribution  of  stress  in  the 
plastic  loacLng,  elastic  unloading  and  plastic  reloading  regions,  as  functions  of  the  strain¬ 
hardening  parameter^ 


1.  Introduction 


Knowledge  of  the  stress  and  deformation  fields  near  the  tip  of  an  advancing  crack  is 
central  to  the  continued  development  of  engineering  fracture  mechanics,  and  important 
advances  have  been  made  in  recent  years.  Under  the  assumption  of  quasi-static  growth  and 
small-strain  ideal  elasto-plasticity,  rigorous  asymptotic  solutions  have  been  found  in  anti¬ 
plane  strain  by  Chitaley  and  McClintock  (1971);  in  plane  strain  by  Slepyan  (1974),  Gao 
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(1980),  Rice  et  al.  (1980)  and  Drugan  etal.  (1982);  and  in  plane-stress  mode  II  by  Ponte 
(1985a).  In  addition,  a  general  discussion  of  the  subject  is  found  in  Rice  (1982). 

The  corresponding  work  for  hardening  elastic-plastic  solids  is  at  a  less  developed 
stage.  Slepyan  (1973)  provided  the  near-tip  fields  of  a  mode  III  crack  propagating  steadily 
in  a  linear-hardening  J2-deformation  theory  material.  Also  for  linear-hardening  materials, 
Amazigo  and  Hutchinson  (1977)  produced  steady-state  asymptotic  solutions  using  J2-flow 
theory  for  modes  1  (plane  strain  and  plane  stress)  and  III.  These  solutions,  however, 
neglected  the  possibility  of  reverse  plastic  flow  on  the  crack  flanks  which  is  known  to 
occur  in  the  elastic  perfectly-plastic  solutions  mentioned  above.  Zhang  et  al.  (1983) 
included  reloading  (as  well  as  a  Bauschinger  effect)  in  their  solution  of  the  plane-strain 
mode  I  problem.  More  recently,  Ponte  (1985b)  provided  a  complete  set  of  solutions 
including  reloading  for  modes  I  (plane  strain  and  plane  stress),  II  (plane  strain  and  plane 
stress)  and  HI. 

Further  work  on  the  hardening  elastic-plastic  steadily  moving  crack  include 
publications  by  Lo  and  Pierce  (1981)  who  considered  a  mode  in  crack  moving  in  a  linear- 
hardening  comer  theory  material,  and  Gao  etal.  (1983)  who  investigated  a  mode  III  crack 
propagating  in  a  power-hardening  J2-flow  theory  material. 

Even  though  deformation  theory  is  unreliable  for  non-proportional  loading  histories 
such  as  that  which  a  near-tip  particle  experiences  as  it  moves  past  a  crack  tip,  it  is  found  by 
comparing  the  solutions  of  Slepyan  (1973),  and  Amazigo  and  Hutchinson  (1977)  for  the 
mode  m  case  that  the  deformation  theory  solution  presents  the  same  general  features  as  the 
flow  theory  solution.  In  fact,  some  of  the  detailed  results  such  as  the  dependence  of  the 
strength  of  the  singularity  on  the  hardening  parameter  are  in  remarkably  close  agreement. 
On  the  other  hand,  for  a  linear-hardening  material,  deformation  theory  suitably  modified  to 
account  for  unloading  and  reloading  yields  an  analytically  much  simpler  set-up  than  flow 
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In  this  work  we  propose  to  solve  the  plane-stress  and  plane-strain  mode  I 
asymptotic  problems  using  J2  deformation  theory  with  linear  strain-hardening,  and  to 

compare  the  results  to  the  corresponding  flow  theory  results.  We  will  make  use  of  concise 
complex-variable  formulations  that  will  yield  straightforward  analytical  solutions  for  the 
near-tip  fields. 

2.  Formulation 

With  reference  to  Figure  2.1,  let  x2  (i  =  1,2,3)  be  a  Cartesian  coordinate  system  of 
fixed  orientation  travelling  with  the  crack  tip  such  that  the  x3  axis  coincides  with  the  straight 
crack  front  Also  let  e(  be  the  unit  vector  corresponding  to  the  xj  direction.  Similarly,  let  r, 
0  be  polar  coordinates  corresponding  to  xa  (a  =  1,2),  and  er  ,  ee  be  the  corresponding 
unit  vectors.  The  crack  tip  moves  with  velocity  V*Vej  with  respect  to  the  stationary 
coordinate  system  X2 .  In  our  steady-state  analysis  the  crack  tip  speed  V  is  constant  so  that 
the  material  derivative  is  given  by 

( )'  -  -V(  ),|  (2.1) 

For  mode  I  loading  of  an  isotropic  material  with  the  given  geometry  we  have  that 
the  only  non-zero  components  of  the  stress  tensor  are  oap  (a,(i  =  1,2)  (and  o33  for  plane 

strain).  Then,  neglecting  inertia,  the  equilibrium  equations  are 

°afca  =  °  (2.2) 

The  only  non-zero  components  of  the  strain  tensor  are  eap  (and  e33  for  plane 

stress).  The  equations  of  compatibility  are 

eaP  ^aX,JJn  =  ®  (2.3) 

where  eap  is  the  two  dimensional  permutation  symbol  (e2  j  =  e^  =  0,  e12  =  -e^  =  1). 

The  constitutive  relations  are 


=  Cy*1  +  EyP* 


(2.4) 


where  the  components  of  the  elastic  strain  tensor  are  given  by 

Eij'i  =  (1/E)  ((1+  v)  a,j-v  o,*  8^]  (2.5) 

and  the  plastic  strain  is  given  by  a  deformation  theory  of  plasticity  modified  to  account  for 
unloading  and  reloading.  From  previous  work  we  expect  to  have  three  different  types  of 
sectors  asymptotically  into  the  crack  tip.  As  depicted  in  Figure  2.2,  we  expect  to  have  a 
plastic  loading  sector  ahead  of  the  crack  (0  <  0  <  0t),  followed  by  an  elastic  unloading 
sector  (0i  <  0  <  02),  followed,  in  some  cases,  by  a  plastic  reloading  sector  (02  <  0  <  rt). 

Thus  when  the  material  particle  P  is  loading  ahead  of  the  crack,  the  components  of 
its  plastic  strain  tensor  are  given  by 

CyP*  *  (3/2)  (1/E,-  1/E)  Sy  (2.6) 

where  Sy  are  the  components  of  the  stress-deviator  tensor,  and  E(  *  Et  (ae)  is  the  secant 
modulus,  which  is  a  function  of  the  effective  stress  oe  =  [(3/2)sijSij]1/2.  For  a  bilinear 
material  with  the  uniaxial  stress-strain  relation  depicted  in  Figure  2.3  equation  (2.6)  reduces 
to 

£ijpi  =  (3/2)  (1/E)  (1/ot-l)  (1-C V0e)  Sij  (2-7) 

where  o0  is  the  yield  stress  in  tension,  and  a  is  the  ratio  of  the  tangent  modulus,  Et,  to  the 
elastic  modulus,  E. 

Note  that  for  |o0/oe|  «  1,  equation  (2.7)  simplifies  to 

EyP1  =  (3/2)  (1/E)  (1/cc-l)  Sy  (2.8) 

Observing  that  for  small  non-dimensional  r,  we  expect  G^a0  =  0(r*)  (where  -1/2  <:  s  <  0 
for  0  <  a  1  1),  it  follows  that,  sufficiently  close  to  the  crack  tip  (0  <  r  «  e’1/l$l),  and 
provided  that  ac  *  0,  equation  (2.8)  is  a  valid  approximation  to  (2.7).  Note  that  the 
approximation  breaks  down  as  a  -»  0.  We  conclude  that  for  0  <  a  £  1 

Ejj  =  (1/E*)  [(1+v*)  Gy  -  v*att  8y]  (2.9) 

where  E*  =  ocE  and  v*  =  a(2v-l)/2  +1/2. 

When  the  particle  unloads  elastically  at  location  1  in  Figure  2.2,  its  plastic  strain  is 
assumed  to  be  continuous  across  the  boundary.  As  it  moves  through  the  elastic  sector,  its 
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plastic  strain  remains  constant  Thus  the  plastic  strain  in  the  elastic  unloading  sector  is 
specified  by 

EijP1  =  [e/]1  (2.10) 

Note  that  [EyP1]1  is  a  function  of  x2  only,  and  is  evaluated  for  a  given  particle  by  letting  xx 
=  x2/tan  0j  in  (2.8). 

When  the  particle  reloads  plastically  at  location  2  in  Figure  2.2,  its  plastic  strain  is 
again  assumed  to  be  continuous  across  the  boundary,  and  hence  is  the  same  as  the  plastic 
strain  that  it  had  just  before  it  unloaded.  Thus  the  plastic  strain  in  the  plastic  reloading 
sector  is  given  by 

EyPi  =  (3/2)  (1/E)  (1/ot-l)  Sy  -  [EyPi]2  +  [CyP1]1  (2.11) 

where  [EyP1]2  is  evaluated  for  a  given  particle  by  letting  =  x2/tan02  in  (2.8) .  Note  that 
any  further  unloading/reloading  can  be  modeled  in  a  similar  manner. 

Equations  (2.2),  (2.3)  and  (2.4)  together  with  (2.5)  and  either  (2.8),  (2.10)  or 
(2.11)  (depending  on  whether  loading,  unloading  or  reloading  occurs,  respectively)  are  the 
governing  equations  in  each  sector. 

We  solve  these  equations  in  any  given  sector  by  first  introducing  a  stress  function  / 
to  satisfy  equilibrium  so  that 

°oP  =  caXePn  (212) 

Next  note  that  in  the  plastic  loading  sector  we  effectively  have  an  elastic  problem 
with  modified  elastic  constants  E*  and  v*,  and  compatibility  yields 

V4/  =  0  (2.13) 

where  we  have  used  the  fact  that  for  plane  stress  o33  =  0,  and  for  plane  strain  e33  =  0 
(which  implies  c33  =  v’o^).  The  solution  of  (2.13)  can  be  expressed  as 

/  =  Re{  z  $(z)  +  x(z)}  (2. 14) 

where  $  and  %  are  analytic  functions  of  z  =  reie  in  the  sector.  Since  we  are  interested  in 
asymptotic  results  we  will  assume  the  following  forms  for  <J>  and  x 


<|>(z)  =  az1 


V(z)  =  X'(z)  =  bz* 


(2.15) 


where  a  and  b  are  undetermined  complex  constants.  The  stress  components  oap  can  then  be 
found  in  terms  of  a,  b  and  s  by  means  of  the  Muskhelishvili  complex  variable  relations 
given  in  Appendix  A. 

In  the  elastic  unloading  sector  compatibility  and  either  the  plane  stress  or  plane 
strain  condition  lead  to 

(l-*2)  V«/  =  -E{  [e„P»]»  4  y  [e33p‘]'  }t22  (2.16) 

where  y  =  0  for  plane  stress,  or  y  =  v  for  plane  strain.  The  solution  to  this  P.D.E.  can  be 
written  as 

/  =  /h  +  /P«  (2.17) 

where  /p*  is  any  particular  solution  of  the  P.D.E,  and  the  homogeneous  solution  /h  is 
given  by  an  expression  similar  to  (2.14)  with  new  functions  $  and  x  (analytic  in  the  sector) 
which  are  of  the  form  of  (2.15)  but  with  different  complex  coefficients  c  and  d.  It  follows 
that  we  can  write 

<*ap  =  <W  +  °aPp“  (2.18) 

where  the  oaph  can  be  expressed  through  /h  in  terms  of  the  new  $  and  y,  and  /p*  can  be 

chosen  such  that 

onP*  -  — [E/(l— y2)]  { [enP!]1  +  y  [e33p1]1  > 

(2.19) 

Cj3p*  =  q77P»  =  0 

In  the  plastic  reloading  sector  compatibility  and  either  the  plane  stress  or  plane  strain 
condition  lead  to 

[1  -  (Y*)2]  V</  =  E*{[  enPi]2  -  [  eup»]1  +  y*  ([e33P']2  -  [e33P«]i)  },22  (2.20) 


This  equation  is  of  the  same  form  as  equation  (2.16),  and  can  be  solved  in  the  same 
manner.  Equation  (2. 18)  applies  with  new  constants  e  and  f  for  the  functions  <!>  and  \|/  in  the 
expressions  for  the  homogeneous  stresses,  and  equations  similar  to  (2.19)  for  the  particular 


stresses. 
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Noting  that  the  inhomogeneous  terms  in  equations  (2.16)  and  (2.20)  are  functions 
of  x2  only,  it  pays  to  consider  the  material  derivatives  of  these  equations  (and  also  of 

equation  (2.13)  for  consistency).  Thus  we  find  that  the  same  homogeneous  equation 
governs  in  all  three  sectors 

V43  =  0  (2.21) 

where  3  =  /'  is  the  stress-rate  function.  The  solution  to  (2.21)  is  given  by 

3  =  Re{  z<D(z)  +  X(z)}  (2.22) 

where  0  and  X  are  analytic  in  each  sector.  For  consistency,  in  the  plastic  loading  sector  we 
let 

O(z)  =  Az*-‘  'F(z)  =  X'(z)  =  Bz*-1  (2.23) 

where  A  and  B  are  undetermined  complex  constants,  but  are  related  to  a  and  b  through 

A  =  -sVa  B  =  -sV(a+b)  (2.24) 

Note  that  even  though  it  will  turn  out  that  it  is  most  convenient  to  use  the  stress-rate 
function  coefficients  A  and  B  in  the  formulation  of  the  problem,  with  this  simple 
conversion  we  can  always  go  back  to  our  original  expressions  above  to  find  ihr  stress 
components.  Expressions  analogous  to  (2.23)  and  (2.24)  apply  in  the  elastic  unloac  g  and 
plastic  reloading  sectors,  but  with  complex  coefficients  pairs  C,  D  and  E,  F,  respectively. 

Later  we  will  need  expressions  for  the  traction-rate  force  F  and  the  velocity  vector  v 
in  terms  of  the  Muskhelishvili  rate  functions,  and  these  are  also  given  in  Appendix  A. 

It  is  clear  from  the  formulation  in  terms  of  the  stress  function  that  this  is  a  fourth 
order  problem.  Given  the  variable-separable  nature  of  the  solution  (z  *  reie),  we  need  to 
impose  four  boundary  conditions  in  8.  Mode  I  symmetry  requires 

u2(r,0)  =  O12(r,0)  =  0  (2.25) 

and  the  traction-free  condition  on  the  crack  faces  requires 

a12(r,Jt)  =  a22(r,tc)  =  0  (2.26) 

Note  that  these  four  conditions  are  homogeneous,  and  lead  to  an  eigenvalue 
problem  for  s  as  a  function  of  a  (and  v).  Thus  we  expect  to  determine  the  stress  and 
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deformation  distributions  only  up  to  an  arbitrary  amplitude  factor  analogous  to  the  stress 
intensity  factor  in  the  elastic  problem. 

Next  we  need  to  connect  the  solutions  in  the  different  sectors  through  some 
appropriate  continuity  conditions.  Thus  we  impose  four  conditions  at  each  unloading  and 
reloading  boundary.  Let  [  ]  denote  the  jump  in  a  quantity  as  0  increases  infinitesimally;  then 
the  traction  components  of  the  stress  must  be  continuous 

[cy  =  [Oee]=0  (2.27) 

and  the  displacement  Held  must  also  be  continuous 

[Uj]  =  [ue]  =  0  (2.28) 

Note  that  in  the  previous  development  we  imposed  continuity  of  the  plastic  strains; 
this  should  be  viewed  as  part  of  the  constitutive  model,  and  not  as  continuity  conditions 
independent  of  (2.27)  and  (2.28). 

The  boundary  and  continuity  conditions  described  above  when  expressed  in  terms 
of  stress  functions  are  fairly  cumbersome  due  to  the  so-called  particular  terms.  It  is  then 
useful  to  rewrite  these  conditions  in  terms  of  rate  quantities.  Thus  the  boundary  conditions 
(2.25)  and  (2.26)  take  the  form 

v2(r,0)  *  X12(r,0)  =  0  (2.29) 

X12(r,Jt)  =  X22(r,n)-0  (2.30) 

where  X  =  O',  and  v  =  u'.  The  last  of  these  two  sets  of  conditions  can  be  rewritten  further 
in  terms  of  the  traction-rate  force 

F^r.Jt)  =  F2(r,Jt)  =  0  (2.31) 

Given  conditions  (2.27),  (2.28)  and  the  continuity  of  the  plastic  strain  across  each 
elastic -plastic  boundary,  it  can  be  shown  (Ponte,  1985b)  that  for  a  strain-hardening  material 
(0  <  a  £  1)  the  stress  tensor  must  be  fully  continuous.  It  then  follows  from  the  equilibrium 

equations  that  the  traction  components  of  the  stress-rate  tensor  must  be  continuous  as  well 
[Xre]  =  [Xee]=0  (2.32) 

This  in  turn  implies  continuity  of  the  traction-rate  force  vector 
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[F1]  =  [F2]  =  0  (2.33) 

Similarly,  it  can  be  shown  (Ponte,  1985b)  that  under  the  given  assumptions  the  velocity 
vector  must  also  be  continuous 

[vj  =  [v2]  =  0  (2.34) 

Thus  conditions  (2.33)  and  (2.34)  will  be  used  instead  of  the  more  primitive  conditions 
(2.27)  and  (2.28)  to  connect  the  fields  in  the  different  sectors. 

Finally,  we  need  to  characterise  the  onsets  of  unloading  and  reloading.  For  flow 
theory,  ae‘  is  continuous  across  the  unloading  boundary,  and  the  unloading  angle,  0!,  is 

uniquely  determined  by  the  condition 

ot'(r,  0!)  =  0  (2.35) 

For  deformation  theory,  however,  oe'  is  not  in  general  continuous.  This  is  due  to  the  fact 
that  the  component  of  the  strain-rate  tensor  Dn  (D  =  E  )  does  not  necessarily  vanish  on  the 

plastic  side  of  the  unloading  boundary  (as  is  the  case  for  flow  theory).  Consequently,  there 
may  be  a  discontinuity  in  1^,  and  therefore  in  <Je‘.  Hence  there  is  a  range  of  possible 

unloading  conditions  with  extreme  conditions 

oe'(r,  Gj-)  =  0  (2.36a) 

Oe(r,  01+)  =  0  (2.36b) 

Of  all  the  possibilities,  condition  (2.36a)  makes  most  sense  physically  since  it  implies  that 
the  effective  plastic  strain-rate  vanishes  at  unloading.  On  the  other  hand,  if  we  think  of 
deformation  theory  as  a  model  for  flow  theory  any  other  unloading  condition  would  be  just 
as  adequate.  Slepyan  (1973)  preferred  condition  (2.36b)  because  it  gave  the  most 

conservative  estimate  for  the  strength  of  the  singularity. 

The  reloading  angle,  02,  is  determined  by  the  condition  that  the  material  particle 


regain  the  effective  stress  level  it  had  when  it  unloaded 
oe(  x2/tan0!,  x2)  =  Oe(  x2/tan  02,  x2) 


(2.37) 
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Note  that  unlike  the  boundary  and  continuity  conditions,  the  unloading  and 
reloading  conditions  (2.36)  and  (2.37)  must  include  particular  terms  since  they  involve 
non-rate  quantities,  and  cannot  be  simplified  any  further. 

Conditions  (2.29),  (2.31),  (2.33),  (2.34),  (2.36a)  and  (2.37)  are  expressed  in 
Appendix  B  in  terms  of  the  eigenvalue  (s),  the  complex  coefficients  of  O  and  (A,  B,  C, 
D,...),  and  the  unloading  and  reloading  angles  (0j  and  02)  for  both  the  two-region  problem 

(no  reloading  occurs)  and  the  three-region  problem  (if  reloading  occurs).  Note  that  for 
plane  stress  these  conditions  are  independent  of  Poisson's  ratio  (v),  but  not  so  for  plane 
strain.  The  solution  procedures  used  to  obtain  the  unknowns  s,  0lf  02,  A,  B,  C,  D,...  are 

also  sketched  in  this  appendix.  Once  these  are  found,  the  stress  distribution  in  the  different 
sectors  is  determined  using  the  expressions  for  the  stresses  determined  previously. 

3.  Results  and  Discussion 

3.1  Plane  stress 

The  results  for  this  case  arc  summarised  in  Figures  3.1  to  3.5  and  in  Table  3.1. 
Figure  3.1  shows  plots  of  the  strength  of  the  singularity,  s,  versus  the  square  root  of  the 
hardening  parameter,  a,  for  two  deformation  theory  materials  satisfying  the  two  extreme 
unloading  conditions  given  by  (2.36a)  and  (2.36b)  in  the  formulation  section,  and  also  for 
the  flow  theory  material  (Ponte,  1985b).  It  is  found,  as  Slepyan  (1973)  observed  for  the 
mode  HI  problem,  that  the  deformation  theory  with  unloading  condition  (2.36b)  results  in  a 
stronger  singularity  (  and  smaller  angular  extent  of  the  plastic  zone)  than  its  counterpart 
with  unloading  condition  (2.36a),  and  that  deformation  theory  in  general  produces  a 
stronger  singularity  than  flow  theory.  It  is  also  seen  that  s  --a1/2  for  small  a. 

Table  3.1  presents  values  of  s,  0j  and  02  for  several  values  of  a  ranging  between 

zero  and  unity  for  the  two  deformation  theories,  and  for  flow  theory.  It  is  seen  that  the 


-11- 


singularities  agree  very  well  for  the  deformation  and  flow  theories,  especially  for  large 
strain-hardening.  For  the  less  singular  deformation  theory  material  (corresponding  to 
condition  (2.36a) )  the  unloading  angle  6j  increases  from  the  elastic  value  to  a  maximum 
value  of  90°  as  a  ->  0.  The  opposite  behaviour  is  observed  for  the  more  singular  material, 
and  ©!  gets  progressively  smaller  as  a  0.  The  flow  theory  material  shows  an 
intermediate  behaviour:  first  increases  as  a  decreases,  and  then  decreases  as  a  -»  0. 

Reloading  occurs  for  the  less  singular  deformation  theory  material  at  the  value  of  the 
hardening  parameter  a*  =  0.215,  it  never  occurs  for  the  more  singular  material,  and  takes 
place  at  a*  -  0.01  for  the  flow  theory  material.  The  actual  size  of  the  reloading  sector  is 
considerable  for  the  less  singular  deformation  theory  material  (up  to  20°),  and  very  small 
for  the  flow  theory  material  (less  than  0.001*). 

Figure  3.2  shows  plots  of  the  particle  (fixed  x2)  effective  stress  for  several  values 

of  a,  and  corresponds  to  the  less  singular  material.  By  comparison  to  the  elastic  result 
plasticity  introduces  two  new  features:  unloading  with  a  discontinuous  slope  [oe  ]  <  0,  and 
the  production  of  residual  stress  on  the  wake  of  the  crack.  These  phenomena  become  more 
significant  for  smaller  a,  and  in  fact  the  need  for  the  reloading  sector  becomes  evident  for 
small  enough  a.  We  also  note  that  this  material  has  two  local  minima  in  the  elastic 
unloading  sector,  whereas  the  flow  theory  material  has  only  one.  Finally,  note  that  for  very 
small  a  the  effective  stress  level  in  the  plastic  zones  approaches  a  constant  value,  thus 
approximating  a  flat  variation  as  expected. 

Figure  3.3  shows  plots  similar  to  Figure  3.2  but  correspond  to  the  more  singular 
material.  It  presents  the  same  general  features  as  the  previous  material  except  for  the  fact 
that  reloading  never  occurs,  and  that  the  behaviour  in  the  elastic  unloading  zone  resembles 
more  closely  that  of  the  flow  theory  material . 

Figure  3.4  shows  plots  of  the  angular  (fixed  r)  variations  of  the  stresses  for  large 
and  moderate  strain-hardening  for  the  deformation  material  with  unloading  condition 
(2.36a).  The  stress  distributions  are  close  to  the  corresponding  flow  theory  distributions 


for  large  strain-hardening,  but  get  progressively  farther  away  from  the  flow  theory 
distributions  as  a  0.  The  main  qualitative  differences  are  the  jumps  in  the  slopes  of  an 
and  oe,  and  the  fact  that  all  the  components  of  the  stress  tensor  approach  0  =  n  with  slopes 
opposite  in  sign  to  their  flow  theory  counterparts.  This  effect  is  most  noticeable  for  an 
which  becomes  positively  unbounded  near  0  =  n,  whereas  it  becomes  negatively 
unbounded  for  flow  theory. 

Figure  3.5  shows  plots  of  the  angular  variations  of  the  stresses  for  large  and 
moderate  strain-hardening  for  the  deformation  material  with  unloading  condition  (2.36b). 
Similar  observations  are  made  about  the  stress  distributions  as  in  Figure  3.4,  but  it  is  noted 
that  in  this  case  the  stress  distributions  agree  more  closely  with  the  flow  theory  results  in 
the  elastic  unloading  sector. 

3.2  Plane  strain 

As  it  was  pointed  out  in  the  formulation  section,  the  results  for  this  case  are 
dependent  on  Poisson's  ratio  (v).  Here  we  chose  not  to  investigate  this  dependence,  and  all 
the  results  presented  in  this  section  are  for  v  =  1/3.  These  results  are  summarised  in  Figures 
3.6  to  3.8,  and  in  Table  3.2. 

Figure  3.6  shows  plots  of  s  versus  a1/2  for  the  deformation  theory  material  with 
unloading  condition  (2.36a),  and  for  the  flow  theory  material.  (For  the  deformation  theory 
material  with  unloading  condition  (2.36b)  no  consistent  solutions  were  found  even  for 
large  strain-hardening).  For  large  strain-hardening  the  two  plots  are  coincident,  but  as  a 
becomes  smaller  the  two  diverge  with  the  deformation  theory  material  showing  the  stronger 
singularity.  For  very  small  a,  the  results  are  unexpected  for  both  materials:  the  deformation 
theory  material  has  a  discontinuous  s  versus  a  curve,  and  the  flow  theory  material  has  a 
non-vanishing  s  as  a  -» 0. 

Table  3.2  presents  values  of  s,  0!  and  02  for  several  values  of  a  for  the  two 
materials  mentioned  above.  It  also  presents  the  values  of  03  and  04  corresponding  to  the 
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beginning  and  end  of  an  isolated  plastic  sector  appearing  inside  the  original  elastic 
unloading  sector  of  the  deformation  theory  material  for  small  values  of  a.  The  deformation 
theory  material  is  found  to  reload  (on  the  crack  flank)  at  a  higher  level  of  hardening  (a*  = 
0.391)  than  the  flow  theory  material  (a*  -  0.3).  Further,  the  three-region  solution  becomes 
inconsistent  for  the  deformation  theory  material  for  0  <  a  <  a+  (a+  =  0.057)  due  to  the 
fact  that  the  material  particle  reloads  inside  the  elastic  unloading  sector,  thus  requiring  a  five 
region  solution  for  0  <  a  <  a+.  The  transition  from  the  three-region  solution  to  the  five- 
region  solution  is  not  smooth  because  the  (discontinuous)  unloading  constraint  produces  a 
redistribution  of  the  stress  and  deformation  fields,  and  it  accounts  for  the  discontinuous  s 
versus  a  curve.  Note  that  appropriate  changes  and  additions  to  the  formulation  presented 

above  are  needed  to  treat  the  five-region  problem,  but  they  are  quite  messy  and  are  not 
presented  here.  The  unloading  angle  0!  is  seen  to  increase  as  a  -4  0  for  both  materials; 

however,  for  deformation  theory  it  has  a  maximum  value  of  90*,  whereas  for  flow  theory  it 
reaches  much  larger  values.  This  fact  might  account  for  the  necessity  of  the  extra  plastic 
loading  sector  for  deformation  theory  as  a  -»  0.  The  reloading  angle  02  is  found  to  be  of 

comparable  magnitude  for  both  cases. 

Figure  3.7  depicts  plots  of  the  particle  effective  stress  for  several  values  of  a. 
Unloading  is  observed,  and  reloading  on  the  crack  flanks  takes  place  at  an  intermediate 
value  of  a.  Note  that  the  unloading  sector  develops  a  bump  which  gets  higher  as  a 
decreases,  and  for  the  critical  value  of  a  -  a+  the  three-region  solution  gives  way  to  the 
five-region  solution.  For  very  small  a,  and  ignoring  the  sharp  unloading  sector,  this  curve 
resembles  the  curve  corresponding  to  the  flow  theory  material. 

Figure  3.8  shows  plots  of  the  angular  variations  of  the  stresses  for  large  and 
moderate  strain-hardening.  We  observe  discontinuities  in  the  slopes  of  ow  O33,  and  oe  at 

both  the  unloading  and  reloading  boundaries.  By  comparison  to  the  stress  distributions  of 
the  flow  theory  material  we  find  that  the  distributions  look  roughly  similar  near  the  crack 


faces,  but  become  quite  different  ahead  of  the  crack  as  a  decreases.  Also  the  deformation 
theory  material  does  not  show  as  high  a  level  of  triaxiality  as  the  flow  theory  material. 

4.  Conclusions 

For  a  linear-hardening  material,  a  deformation  theory  of  plasticity,  adequately 
modified  to  take  into  account  unloading  and  reloading,  yields  an  approximate  closed  form 
solution  to  the  stress  and  deformation  fields  near  the  tip  of  a  steadily  propagating  crack. 
Inside  a  region  of  size  r  =  e*1/W,  the  constitutive  relations  can  be  linearised  to  yield 
modified  linear  elasticity  problems  in  the  different  sectors  that  can  be  solved  analytically 
using  concise  complex-variable  methods.  The  range  of  validity  of  the  approximation 
vanishes  as  the  hardening  parameter  (and  therefore  the  singularity)  becomes  smaller. 

For  large  strain-hardening,  the  functional  dependence  of  the  singularity,  s,  on  the 
hardening  parameter,  a,  for  deformation  theory  is  very  close  to  that  for  flow  theory  in  both 
plane  stress  and  plane  strain.  Remarkably,  the  s  versus  a  curves  for  deformation  and  flow 
theory  remain  close  for  intermediate  and  small  strain-hardening  for  the  plane  stress  case 
(and  the  anti-plane  shear  case  as  pointed  out  by  Amazigo  and  Hutchinson,  1977),  but  not 
for  the  plane  strain  case  where  unexpected  results  are  obtained  for  both  theories.  In  general, 
however,  deformation  theory  provides  a  stronger  singularity  than  flow  theory  in  all  cases 
studied  (including  mode  III). 

As  expected,  for  large  strain-hardening  the  angular  variations  of  the  stresses  are 
similar  for  both  theories,  but  they  become  quite  different  as  the  hardening  becomes  smaller. 
Thus  the  behaviour  of  the  sizes  of  the  plastic  and  elastic  sectors  as  measured  by  and  02 
are  also  fairly  different  fra1  both  theories  as  a  — >  0. 

While  for  flow  theory  unloading  is  uniquely  defined  by  the  condition  oe'(r,0,)  =  0, 
for  deformation  theory  we  must  admit  a  discontinuity  in  <?e'  at  unloading  boundaries,  and 


there  is  a  range  of  possible  unloading  conditions  between  the  two  extreme  conditions 
oe’(r,0i-)  =  0  and  ce'(r,0i+)  =  0.  It  was  found  that  the  first  condition  leads  to  a  plastic 

loading  zone  ahead  of  the  crack  of  larger  angular  extent  than  that  of  the  second,  and  also  to 
a  weaker  singularity  s  for  a  given  value  of  a. 
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Appendix  A.  Muskhelishvili  formulation 


Stress 


Given  the  Muskhelishvili  functions  <j>  and  \y,  we  can  express  the  stress  components 
corresponding  to  a  homogeneous  elastic  problem  in  the  following  fashion 
ctn  +  c22  =  2[  <>'(  z)  +  $’(  z)J 


(A.l) 


cr22-  On  +  2ic12  =  2[  z  *"<z)  +  v'(z)] 


Traction-rate  force 

The  traction-rate  force  on  a  boundary  C  with  fixed  origin  z„  and  variable  endpoint  z 
is  given  by 

F1  +  iF2  =  Jc  [T^s)  +  iT2(s)]  ds  (A.2) 

where  s  is  the  arc  length  measure  along  C,  Ta=  oapnp  are  the  traction  stress-rates  and  np 

define  the  unit  normal  to  C. 

The  integral  in  (A.2)  is  path-independent,  and  it  can  be  shown  that 

F1  +  iF2  =  -i(3§1  +  i3,2)  (A.3) 


Velocity 

The  velocity  expressions  depend  on  the  sector 

(i)  Elastic  unloading  sector 

2G(vt  -  iv2)  =  K  <D(  z)  —  z  <D’(z)  -  'F(z)  (A.4a) 

where  k  =  (3-v)/(l+v)  for  plane  stress,  or  k  =  3-4v  for  plane  strain. 

(ii)  Plastic  loading  sectors 

2G*(vt  -  iv2) .  k*  d>(  z)  -  z  <D'(z)  -  'P(z) 
where  k*  =  (3-v*)/(l+v*)  for  plane  stress,  or  k*=  3-4v*  for  plane  strain. 


(A.4b) 
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Appendix  B.  The  governing  conditions  in  terms  of  the  unknown  constants 
Two-region  problem  (a*  <  a  <  1). 

In  this  problem  we  have  six  unknowns:  two  real  (s  and  ©j),  and  four  complex  (A, 
B,  C,  D)  for  a  total  of  ten  real  unknowns.  The  boundary  and  continuity  conditions  reduce 


Symmetry  (0  =  0) 
A  =  A 


B  =  B 


Vanishing  traction-rates  (8  =  n) 
s  C  +  ei2ntC  +  D  =  0 
Continuity  of  traction-rates  (0  =  0^ 

sea0i  (  C-  A)  +ei2*9l  (C -  A)  +  (  D-  B)  =  0 
Continuity  of  velocities  (8  =  0j) 

(1+v*)  (1+K*)  A  -  a(l+v)  (1+k)  C  + 

(3/2)  (a-1)  [A  +  sei2(*'1)eiA  +  e^iB]  =  0 
The  unloading  condition  (oe  (r,0i-)  =  0)  implies 
[(4/3)  Y*  COS(S0j)  008(5-1)0!]  A2  + 


(B.l) 


(B.2) 


(B.3) 


(B-4) 


(B.5) 


cos  ©1{(s-l)  [s+3-4cos20i]  A2  +  (s-1)  (4cos20i~2)  BA  +  B2}  =  0 
where  y*  =  1  for  plane  stress,  or  y*  *  1  +  4v*(v*-l)  for  plane  strain. 

Given  s,  0j  and  (B.l),  equations  (B.2)  to  (B.4)  can  be  viewed  as  a  system  of  six 
linear  homogeneous  equations  in  the  six  unknowns  A,  B,  C,  C,  D  and  D.  The  condition 

for  a  non-trivial  solution  is  that  the  determinant  of  the  system  matrix  be  zero,  or 
[(l-cos20j)  sin20j]  s3  +  [(1-6)  (l-cos20i)  sin(2s0j)]  s2  + 

{(l/2)co(8-co)  +  (1+©)  [cos(2s0i)-l]}  sin(20!)  s  + 


(B.6) 


{(1-8)  (1+©)  [cos2s(01-Jt)-l]  sin(2s0!)  -  (1/2)  ©  (6+©)  sin(2s0j)  + 
(l/2)©8sin[2s(20rn)]  }*  0 

where  ©  =  (2/3)  a  (1+v)  (1+k)  /(l-o)  and  5  =  (2/3)  (1+v*)  (1+k*) /(I -a). 


WAN 


We  can  next  use  any  five  of  the  six  linear  equations  to  solve  for  B,  C  and  D  in 
terms  of  A.  Given  these  we  can  look  at  equation  (B.5)  as  a  non-linear  relation  between  s 
and  6^  which  together  with  equation  (B.6)  provides  a  system  of  two  non-linear  algebraic 
equations  in  s  and  0j  to  be  solved  numerically  for  a  given  a  (and  v).  It  is  observed, 
however,  that  equations  (B.  1 )  to  (B.6)  are  independent  of  v  for  the  plane  stress  case. 
Hence,  for  plane  stress  s,  0j  and  the  stress  distribution  are  all  independent  of  v. 

Three-region  problem  (a?  <  a  <  a*). 

Here  we  have  nine  unknowns:  three  real  (s,  0lf  02  )  and  six  complex  (A,  B,  C,  D, 
£,  F)  for  a  total  of  fifteen  real  unknowns.  The  boundary  and  continuity  conditions  include 
equations  (B.l)  to  (B.4)  plus  two  more  equations  of  the  type  of  (B.3)  and  (B.4)  applied  to 
the  reloading  boundary  (0  =  62).  The  same  unloading  condition  (B.5)  applies,  but  we  must 

also  include  a  reloading  condition  which  for  brevity  will  not  be  stated  here.  In  principle, 
however,  given  s,  0,,  02  and  (B.l),  wc  can  identify  a  system  of  ten  homogeneous  linear 
equations  in  the  ten  unknowns  (B,  B,  C,  C,  D,  D,  E,  E,  F,  F).  The  zero  determinant 
condition,  and  the  unloading  and  reloading  conditions  then  provide  three  non-linear 
algebraic  equations  in  s,  0j  and  02  which  can  be  solved  numerically  for  given  a  (and  v). 
Again  it  is  noted  that  for  plane  stress  s,  0t  and  the  stress  distributions  are  independent  of  v. 


Tables 


Table  3.1  -  Plane  stress  mode  I 


a  s(l) 

s(2) 

s(3) 

Qjd) 

6j(2) 

ej(3) 

e2(i) 

02(2) 

e2(3) 

1.  -0.5 

-0.5 

-0.5 

79.92 

79.92 

79.92 

. 

- 

0.75  -0.469 

-0.469 

-0.468 

81.11 

73.10 

80.78 

- 

- 

- 

0.667  -0.456 

-0.456 

-0.455 

81.58 

70.52 

80.90 

- 

- 

- 

0.5  -0.422 

-0.423 

-0.420 

82.63 

64.92 

80.76 

- 

- 

- 

0.25  -0.340 

- 

-0.335 

84.52 

- 

78.55 

- 

- 

- 

0.1  -0.241 

-0.255 

-0.237 

86.10 

46.75 

73.65 

178.7 

- 

- 

0.01  -0.0896 

-0.107 

-0.0863 

88.28 

32.94 

61.09 

168.8 

- 

180.0 

0.001  -0.0310 

-0.0423 

-0.0287 

89.29 

22.61 

53.20 

160.9 

- 

180.0 

0.0001  -0.0102 

- 

-0.0093 

89.74 

• 

49.00 

156.9 

180.0 

(1)  Deformation  theory  with  unloading  condition  oe‘(r,  0]-)  =  0 

(2)  Deformation  theory  with  unloading  condition  oe‘(r,  0]+)  =  0 

(3)  Flow  theory  with  unloading  condition  oe'(r,  0j)  =  0 


-0.5 

-0.5 

88.62 

88.62 

- 

- 

- 

- 

75 

-0.478 

-0.478 

89.19 

91.91 

- 

- 

- 

- 

67 

-0.468 

-0.469 

89.34 

93.65 

- 

- 

- 

- 

.5 

-0.441 

-0.442 

89.61 

98.40 

- 

- 

- 

- 

,2 

-0.343 

-0.300 

89.94 

113.88 

174.8 

179.7 

- 

- 

1 

-0.274 

-0.197 

89.99 

123.35 

166.0 

173.6 

- 

- 

05 

-0.269 

-0.142 

89.99 

130.61 

139.3 

161.0 

103.78 

137.38 

04 

-0.264 

-0.127* 

90.00 

131.76* 

158.2 

157.1* 

100.93 

141.60 

01 

• 

-0.0797 

• 

135.19 

. 

145.5 

_ 

(1)  Deformation  theory  with  unloading  condition  ac'(r,  0i~)  =  0 

(2)  Flow  theory  with  unloading  condition  ae'(r,  0!)  =  0 
(*)  Interpolated  result 


Figure  3.1 — Strength  of  the  regularity  in  plenc  iimi  for  (i)  deformation  theory  with 
unloading  condition  (2  36a),  (b)  deformation  theory  with  unloading  condition  (2.36b).  and 


(c)  flow  theory. 


Angle.  6 

Figure  3-J» — Angular  (treat  diitribatioa  m  plant  ww  far  deformation  theory  with 
anloading  condition  (2.36b)  for  huge  drain  hardening,  aornwliarri  wch  that  o,(r,8,) »  1. 
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Fijui*  3.6 — Savagth  of  tb*  tiDgiilariqr  ■  piae  Mraia  (v  •  U3)  (or  (a)  deformation  theory 
with  unloadirtf  oandibao  (2J6e),  (b)  flow  theory. 
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